Families of Vortex Solitons in Periodic Media 
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Various families of charge-one vortex solitons in two-dimensional periodic media are reported. 
These vortices reside either in the semi-infinite gap or higher band gaps of the media. For both Kerr 
and saturable nonlinearities (either focusing or defocusing), infinite vortex families are found. All 
these families do not bifurcate from Bloch bands; rather, they turn around before reaching edges 
of Bloch bands. It is further revealed that vortices with drastically different topological shapes can 
belong to the same vortex family, which is quite surprising. 
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I. INTRODUCTION 



Nonlinear wave propagation in periodic media is an 
important phenomena, and it arises in a wide array of 
physical settings such as photonic-crystal fibers [J, ||, 
photorefractive crystals imprinted with a photonic lat- 
tice [1, H|, and Bose-Einstein condensates loaded in an 
optical lattice [1, d, 0). The interplay between nonlin- 
earity and periodicity results in a lot of new physical 
effects of the wave propagation. For instance, in a pe- 
riodic media, solitons can exist not only under focus- 
ingnonlinearity, but also under defocusing nonlinearity 
[ITS E3] ■ A distinctive feature of the periodic media is 
that it can support a wide variety of solitons residing in 
different band gaps. Examples include fundamental soli- 
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1^,JT3|, dipole solitons 0E3], 
2A [23| . reduced-symmetry 
higher-band vortex solitons 0, [26|, l27f . 

_2c| , and so on — many of which 
have been experimentally observed. Solitons in Bessel- 
ring lattices have been reported too [29|, [3(| ■ In an ef- 
fort to classify solitons in two-dimensional (2D) periodic 
media, Shi and Yang [3l| studied low-amplitude solitons 
near Bloch-band edges. These solitons are Bloch waves 
modulated by slowly-varying envelope functions. Using 
asymptotic techniques, they derived envelope equations, 
based on which they successfully classified all soliton fam- 
ilies bifurcating from Bloch-band edges. Examples in- 
clude fundamental solitons, reduced-symmetry solitons, 
dipole-array solitons, and many others. However, a pe- 
culiar question arose on vortex solitons. According to 
the envelope equations, those vortex solitons should also 
bifurcate from Bloch-band edges. However, numerical 
results suggest that they disappear before reaching band 
edges [3XT33j . This paradox has not been resolved yet, 
and it clouds over the vortex-soliton research in periodic 
media. 

In this paper, we comprehensively investigate charge- 
one vortex solitons residing in various band gaps of a 2D 
periodic media, and resolve the above paradox on this 
subject. For both the Kerr (cubic) and saturable nonlin- 
earities, we find that there exist infinite families of on-site 
and off-site vortex solitons in the semi-infinite gap (for 
focusing nonlinearity) and in the first gap (for defocusing 



nonlinearity) . We further show that all these vortex fam- 
ilies do not bifurcate from edges of Bloch bands. Rather, 
before reaching band edges, they turn back and move 
into band gaps again. Within each vortex family, shapes 
of vortex solitons undergo drastic changes as the prop- 
agation constant varies. As a result, vortices with very 
different topological structures can actually belong to the 
same family — a phenomenon which has never been re- 
ported before. These findings significantly deepen our 
understanding on vortex solitons in general periodic me- 
dia. 



II. VORTEX SOLITONS UNDER KERR 
NONLINEARITY 

We first consider vortex solitons in 2D periodic media 
under Kerr nonlinearity. The mathematical model for 
this situation is 

iU z + U xx + U yy + V{x,y)U + a\U\ 2 U = 0, (1) 

where U is a complex function, a takes the value of 1 for 
focusing nonlinearity and —1 for defocusing nonlinearity, 
and V(x, y) is a 2D periodic potential. This equation 
arises naturally in Bose-Einstein condensates loaded in 
an optical lattice [fa, [y, 0] as well as nonlinear light propa- 
gation in laser-written waveguides [34[ and photonic crys- 
tal fibers with weak transverse index variation. Without 
loss of generality, we take the periodic potential V(x, y) 
as 

V(x,y) = -V {sm 2 x + sm 2 y), (2) 

and take Vq = 6 in all our calculations. 

Vortex solitons in Eq. (J]) are sought in the form 

U(x,y;z) = u(x,y)exp(-ifj,z), (3) 

where /i is the propagation constant, and u(x, y) is 
a complex-valued localized function which satisfies the 
equation 



Uxx + u, y y + [u + V(x, y)]u + a\u\ 2 u = 0. 



(4) 



We will determine these solutions numerically using 
the modified squared-operator iteration method and the 
power-conserving squared-operator method developed in 
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A. The Case of Focusing Kerr Nonlinearity 



First, we study vortex solitons supported by the fo- 
cusing Kerr nonlinearity, i.e. a = 1 in Eq. (j4]). Two 
types of vortex solitons — on-site and off-site ones, will 
be sought. On-site vortices are centered on a lattice site 
(lattice peak), while off-site vortices are centered between 
lattice sites (lattice minimum) [l3|, [2(| [2l|, • For each 
type of vortices, various families are found in the semi- 
infinite gap. The power diagrams of the first two fam- 
ilies of on-site and off-site vortex solitons are shown in 
Fig. [T] (a, b) respectively. Here the power is defined 
by P = \u\ 2 dxdy. A distinctive feature about 

these power curves is that each curve has a slanted U- 
shape. At each propagation constant /i away from the 
band edge, one can find two vortex solitons with differ- 
ent powers within each family. Another feature about 
these curves is that none of them reaches the edge of 
the (first) Bloch band. This means that none of these 
vortex families bifurcates from Bloch bands. These fea- 
tures are in stark contrast with those for soliton families 
reported in [3lJ, which do bifurcate from Bloch bands 
(in the low-amplitude limit), and each propagation con- 
stant only corresponds to a single soliton solution within 
each family. To find out what types of vortex solitons 
are contained in each family, we display them at four 
representative locations of each power curve in Fig. 2 
(on-site) and Fig. 3 (off-site). First we examine the first 
on-site vortex family whose power curve is the solid (blue) 
one in Fig. 1(a). At four marked positions 'a-d' on this 
curve, the corresponding vortex profiles are displayed in 
Fig. 2 (a-d) respectively. The phase structures of all these 
vortices are similar to that of the ring vortex soliton in 
bulk media and shown in the inset of Fig. 2(a). Winding 
around the center of the vortex, the phase increases by 
2ir, thus these vortices have charge one. The intensity 
distributions of these vortices, on the other hand, dif- 
fer significantly at various locations of the power curve. 
On the lower branch of the power curve, when /i is far 
away from the Bloch band, the vortex soliton consists of 
four main humps residing in the four lattice sites closest 
to the center site and forming a diamond configuration 
[see Fig. [2] (a)]. Hence this is the familiar on-site vor- 
tex soliton which has been reported before [2(| H3, Hi| . 
When /i increases toward the band edge along the lower 
branch, the vortex's power as well as its peak intensity 
decreases, and the vortex starts to develop tails on the 
outside of its four main humps [see Fig. H](b)]. As \x fur- 
ther increases beyond a cut-off value ^ =4. 107, which is 
close to but below the band edge ^t e( j se =4.126, vortices in 
this family can not be found anymore. This phenomenon 
was first reported by one of the authors in [33], and it 
was puzzling at the time since we did not expect a so- 
lution branch terminating abruptly away from edges of 
Bloch bands. This puzzle is now resolved by Fig. 1(a). 
What happens is that when [i reaches the cut-off value 
fjtc , the power curve turns around and enters the upper 



branch of the same family [see Fig. 1(a)]. When pass- 
ing through the turn-around, the vortex undergoes rapid 
shape changes and evolves into eight main humps in a 
square configuration surrounding the central lattice site. 
Tails at the outside of this square are significant too. 
As /i moves away from the Bloch band along the up- 
per branch, the tails gradually fade away, and the eight 
main humps intensity totally dominate the vortex [see 
Fig. [2] (d)]. Thus, in this first family of on-site vortex 
solitons, the vortex undergoes dramatic shape changes 
from a four-humped diamond- like structure into an eight- 
humped square structure. This is quite unexpected and 
was not known before. 

Regarding the second on-site family, vortex solitons at 
four representative positions 'e-h' of the power curve are 
displayed in Fig. O (e-h) . The phase structures of these 
vortices are similar to those in the first family [see inset 
of Fig. 2(e)], i.e. these vortices also have charge one. 
But the intensity patterns in this second family are quite 
different from those in the first family. Specifically, in 
the lower branch of the second family, away from the 
Bloch band, the vortex contains eight main humps in a 
larger-diamond configuration [see Fig. 2(e)]. As /i moves 
continuously from the lower branch to the upper one, 
the vortex continuously evolves from this eight-humped 
diamond configuration to a sixteen-humped configuration 
resembling the surrounded wall of a castle [see Fig. 2(e- 
h)]. This dramatic shape change within the same solution 
family occurs for all vortex families, and it is one of the 
main findings of this article. 

For off-site vortex solitons in the semi-infinite gap, var- 
ious solution families arc found as well. In the first off-site 
family, representative solution profiles are displayed in 
Fig. 3 (a-d). These vortices also have a simple 2it phase 
winding structure around the vortex center, thus having 
charge one [see inset in Fig. 3(a)]. Regarding their in- 
tensity patterns, when \i is far away from the Bloch band 
on the lower branch, the vortex has four main humps oc- 
cupying four adjacent lattice sites in a compact square 
configuration. This is the familiar off-site vortex soliton 
reported in [l3|, [2l|, [22|, [HJ • As /i moves from the lower 
branch to the upper one, the vortex changes from the 
four-humped square configuration to a twelve-humped 
compact cross configuration [see Fig. 3(a-d)]. In the 
second off-site family, lower-branch vortex solitons away 
from the Bloch band have eight main humps arranged in 
an octagon configuration [see Fig. 3(e)]. As \i moves to 
the upper branch, however, the vortex becomes a square 
with twelve main humps on its perimeter [see Fig. 3(h)]. 
Notice from Fig. [T] (b) that the upper branches of these 
two off-site families have almost the same power for the 
same propagation constant. This is because on these two 
upper branches, vortices of both families have the same 
number of main humps, and nearby humps are all aligned 
along the lattice direction and separated by one lattice 
spacing [see Fig. 3(d,h)]. One may recall that for on- 
site vortices, the upper branch of the first family and 
the lower branch of the second family also have the same 
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FIG. 1: (Color online) Power diagrams of the first two families of on-site (a) and off-site (b) vortex solitons in the semi-infinite 
gap under focusing Kerr nonlinearity. The inset in (a) shows the 2D square lattice of Eq. ((2J); the inset in (b) zooms in on the 
graph near the band edge. Vortex profiles at the marked points are shown in Fig. [2] (on-site) and Fig. [3] (off-site) respectively. 
Shaded: the first Bloch band. 
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FIG. 2: (Color online) Vortex profiles (|lt|) in the first (top) and second (bottom) on-site solution families under focusing Kerr 
nonlinearity at marked points in Fig. [l](a). The inset in (a) is the typical phase structure of all these vortices. The background 
circles represent the lattice sites (with high V values), as in Figs. 3, \5\ and [6] as well. 



number of humps [see Fig. 2(d, e)]. But adjacent humps 
in the latter case are aligned along diagonal directions 
of the lattice and have larger spatial separations than in 
the former case. This caused noticeable power differences 
between these two branches as can be seen in Fig. Q] (a) . 

From Fig. 1, we see that both on-site and off-site vor- 
tices can exist close to the Bloch band [see points 'c, 
g' in Fig. 1(a) for example]. We find that such vor- 
tices near the Bloch band have lower amplitudes. For in- 



stance, the vortex at point 'c' in Fig. 1(a) has amplitude 
0.55, and the one at point 'g' has amplitude 0.48. Com- 
paratively, vortex amplitudes far away from the Bloch 
band are much higher. Additionally, vortices close to 
the Bloch band have long tails, and they resemble the 
Bloch wave of the band edge modulated by a ring- vortex 
envelope [see Fig. 2(c, g) for instance]. This seems to 
suggest that vortex solitons near Bloch bands can be 
treated analytically by asymptotic methods. However, 
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FIG. 3: (Color online) Vortex profiles in the first (top) and second (bottom) off-site solution families under focusing Kerr 
nonlinearity at marked points in Fig. [1] (b) . The inset in (a) is the typical phase structure of all these vortices. 



such asymptotic analysis encounters subtle problems for 
vortex solitons, as we will explain below. Low-amplitude 
Bloch-wave packets in Eq. ([1]) have been analyzed in 31 1 . 
Near the band edge Hedge in Fig. 1 and under focus- 
ing nonlinearity, the leading-order asymptotic solution is 
u(x,y) = eA(X,Y)p(x,y), fi = fi e d 9 e - e 2 - Here, p(x,y) 
is the Bloch wave at the band edge [ledge (with T-point 
symmetry), X = ex, Y — ey, e <C 1, and the envelope 
function A(X, Y) satisfies the equation 



Di(A X x + A YY ) - A + a \A\ 2 A = 0, 



(•5) 



where D\ > is the second-order dispersion coefficient 
at the band edge, and ctQ > is a constant. The envelope 
equation © admits a ring- vortex solution A — /(i?)e ,e , 
where (R, 0) is the polar coordinates of the (X, Y) 
plane. Due to some additional constraints on the solution 
u(x, y), this ring envelope can only be centered at certain 
locations at or between lattice sites HH- Centering the 
ring-vortex envelope at a lattice site, the corresponding 
solution u{x,y) would then be a low-amplitude on-site 
vortex soliton which resembles the ones such as Fig. 2 
(c, g). Centering the ring- vortex envelope between a lat- 
tice site, we would get a low-amplitude off-site vortex 
soliton which resembles the ones in Fig. 3 (c, f). How- 
ever, a contradiction between this asymptotic analysis 
and numerical results is that, according to the asymp- 
totic analysis, these vortices should exist continuously 
as fi approaches the band edge ix e dge (i-e. e — > 0), but 
the numerical results in Fig. 1 indicate that true vortex 
solitons actually do not approach the band edge for each 
vortex family. Apparently, the asymptotic analysis above 
is not entirely correct for vortex solitons, and it must be 



revised in order to explain the true vortex behaviors in 
Figs. 1-3. How this can be done is still unclear yet. Note 
that these vortex behaviors are similar to multi-packet 
solitary waves in the fifth-order Korteweg de Vries equa- 
tion [37]], so the analytical treatment of [37| might be 
useful for the present problem. This remains to be seen. 

Fig- [Hjust shows the first two families of on-site and 
off-site vortex solitons. We have also found other vor- 
tex families in Eq. ^ whose power curves have similar 
slanted U-shapes and are higher than those in Fig. [T] 
Vortices in these higher families contain more intensity 
humps which are located further away from the center of 
the vortex, but their phase structures remain similar to 
those of the first two families (i.e. they all have charge 
one). It can be inferred that infinite families of such vor- 
tex solitons exist in the semi-infinite gap, and all of them 
do not bifurcate from the edge of the first Bloch band. 
Within each higher family, vortex shapes also undergo 
drastic changes as /i moves from the lower branch to the 
upper one. 



The vortex solitons studied above are symmetric about 
the vortex center. Asymmetric vortices also exist in a 
square lattice [23]. We find that the power curves of 
asymmetric vortices also have slanted U-shapes as in Fig. 
1, and they do not bifurcate from Bloch bands either. 
We also find that asymmetric vortices reported in [23| 
belong to the lower branches of these power curves. As 
H moves to the upper branches, shapes of asymmetric 
vortices become very different and more complex. 
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FIG. 4: Power diagrams of the first family of on-site (a) and off-site (b) vortex solitons in the first band gap under defocusing 
Kerr nonlinearity. The insets zoom in on the graphs near the band edge. Vortex profiles at the marked points are shown in 
Fig. [5] (on-site) and Fig. [6] (off-site) respectively. Shaded: the first two Bloch bands. 
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FIG. 5: (Color online) On-site gap vortex profiles (|u|, top row) and phase structures (bottom row) corresponding to the marked 
points in Fig. [4] (a) under defocusing Kerr nonlinearity. Bottom right: color bar of the phase figures. 



B. The Case of Defocusing Kerr Nonlinearity 

In this subsection, we examine vortex solitons under 
defocusing Kerr nonlinearity, i.e. a = — 1 in Eq. (|4|). 
These vortices exist as infinite families in the first band 
gap. Fig. |4]shows the power curves of the first on-site and 
off-site gap-vortex families. It is seen that these power 
curves also have a U-shape, but slanted in an opposite 
direction from the focusing case (see Fig. 1). Spatial 
profiles and phase structures of these gap vortices at the 
marked points in Fig. [5] are displayed in Fig. [5] (on-site) 
and Fig. [5] (off-site). We see that the shapes of these first- 
family defocusing gap vortices (|u|) resemble those of the 



first-family focusing vortices. Specifically, for the first on- 
site defocusing family, as /i moves from the lower branch 
to the upper one, the shape of gap vortices changes from a 
four-humped diamond configuration to an eight-humped 
square configuration (see Fig. [5]), analogous to the fo- 
cusing case (see Fig. 2). For the first off-site defocusing 
family, as [i moves from the lower branch to the upper 
one, gap vortices change from a four-humped square con- 
figuration to a twelve-humped cross configuration (see 
Fig. [6|, also analogous to the focusing case (see Fig. 3). 
The most important difference between defocusing vor- 
tices and focusing vortices lies in their phase structures. 
The phase fields of focusing vortices have a simple 27r- 
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FIG. 6: (Color online) Off-site gap vortex soliton profiles (|u|, top row) and phase structures (bottom row) corresponding to 
the marked points in Fig. |4](b) under defocusing Kerr nonlinearity. Bottom right: color bar of the phase figures. 



winding structure like in ring vortices of bulk media [see 
insets in Figs. 2(a) and 3(a)]. However, phase structures 
of defocusing gap vortices are much more complex (sec 
Figs. [5] and[H]). In the central region of gap vortices, 
the phase field has a simple 27r-winding structure around 
the vortex center. But in the outer region, vortex phases 
at adjacent lattice sites sometimes have it difference [see 
Figs. [Hb,c) and[5]Jb,c)], but some other times do not [see 
Figs, [^d) and|5Ja)]. These complicated phase structures 
make the concept of topological charges ambiguous and 
not well defined for gap vortices. Note that gap vor- 
tices reported earlier in [25|, [33| correspond to the lower 
branches of our gap vortices in Fig. 3J 

From Fig. [?J we see that gap vortices can exist quite 
close to the right edge of the first Bloch band. At this 
band edge, the second-order dispersion coefficient Di 
is negative. Thus according to the envelope equation 
for this band edge, which is similar to Eq. ([5]) with 
only the sign of the A term reversed, ring-vortex so- 
lutions A = f(R)e lB exist under defocusing nonlinear- 
ity (where «o < 0). Hence for fi close to this band 
edge, one may construct the analytical gap-vortex so- 
lution u(x,y) — eA(X,Y)p(x 7 y), where p(x,y) is the 
Bloch wave at this right edge (with M-point symmetry) . 
This analytical gap-vortex solution qualitatively repro- 
duces the intensity distribution, and more importantly, 
the intricate phase structure, of the true solution near 
the right edge of the first band [see Figs. [5jb) and[5Jb) 
for instance]. In particular, it explains the simple 27r 
winding of the phase near the vortex center, and the 7r 
phase difference between lattice sites in the far field (due 
to the M-point symmetry of the Bloch wave at this band 
edge). Of course, this analytical gap- vortex solution can 



not explain why true gap vortices exist as infinite fam- 
ilies, and each family does not continuously reach the 
band edge — the same difficulty we have seen with the 
focusing nonlinearity in the previous subsection. 



III. VORTEX SOLITONS UNDER SATURABLE 
NONLINEARITY 

In photorefractive crystals where many lattice-soliton 
experiments were performed in recent years, the nonlin- 
earity is saturable [38]. Even though wave phenomena 
under saturable and Kerr nonlinearities have many fea- 
tures in common, significant differences exist as well. In 
this section, we examine vortex solitons under saturable 
nonlinearity. 

A simple theoretical model for solitons in photorefrac- 
tive crystals imprinted with a photonic lattice is [H, [38| 



E 



^xx \ U'yy 1,7-/ \ , i 1 9 

l + I{x,y) + \u\ 2 
where E is the applied dc field, 

I(x,y) = I cos 2 (x) cos 2 (y) 



U = —/JiU, 



(6) 



(7) 



is a periodic square-lattice function, Iq is the peak inten- 
sity of this lattice, and [i is the propagation constant. 
Here all variables have been non-dimensionalized [3?j ]. 
When the soliton amplitude is low (|u| <C l), the sat- 
urable nonlinearity becomes the same as the Kerr nonlin- 
earity. Therefore, we can expect that near Bloch bands, 
the vortex behavior under saturable nonlinearity should 
be similar to that under Kerr nonlinearity. In partic- 
ular, they can not bifurcate from Bloch bands either. 
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Below we numerically investigate vortex solitons in the 
semi-infinite gap in Eq. ^ under focusing nonlinearlty 
(Eq > 0). Without loss of generality, we choose param- 
eters Jo = 3 and Eq — 6. Vortex solitons in the first 
band gap under defocusing saturable nonlincarity will be 
briefly described in the end of this section. 

First, we examine on-site vortex solitons under focus- 
ing saturable nonlinearity. An infinite number of such 
vortex families exist in the semi-infinite gap, and one of 
them is reported below. The power curve of this family is 
displayed in Fig. [JJ We see that this curve is quite com- 
plicated. It has a number of turning points (twists), and 
winds up and down in an unexpected way (a somewhat 
similar power curve was reported for ID lattice soltons 
under cubic-quintic nonlinearities in [39]). All vortices 
in this family have simple 2tt phase winding structures, 
thus having charge one [see inset of Fig. Ufa)]. Their 
intensity structures along the power curve are displayed 
in Fig. [5] It is seen that far from the band edge (point 
'a' in Fig. [7]), the power is very high, and the vortex has 
a ring profile [see Fig. HJa)] . This ring profile can be 
easily understood since the vortex here has high inten- 
sities and dominates over the lattice, thus it effectively 
becomes a ring vortex in a lattice-free (bulk) media. As 
the power and intensity decrease (when fi increases) , the 
lattice effect starts to appear, and the vortex reshapes 
itself by distributing its main power to the lattice sites. 
This leads to a four-humped on-site vortex as seen in Fig. 
Etb, c). This vortex has been theoretically analyzed in 
|3q ] and experimentally observed in [22| . This reshaping 
process, however, is a delicate one as is evidenced by the 
two turning points near 'b' in the power curve of Fig. [JJ 
(a). As the power curve gets close to the first band (see 
points 'c, d'), the vortex starts to spread out, and its am- 
plitude becomes low. Vortices here closely resemble those 
of the first on-site Kerr family near the first band [see 
Fig. (b, c)]. Therefore, the power curve of saturable 
vortices can not reach the first band either, and has to 
turn around before reaching the band edge. After the 
turning point, the power starts to increase (as fi moves 
away from the Bloch band), and the vortex reshapes it- 
self again. This time, after a complex reshaping process, 
the eight-humped square vortex at point 'd' becomes a 
four-humped square vortex with its humps residing near 
the four diagonal lattice sites around the vortex center 
[see Fig. [5] (e, f)]. Strangely, in this reshaping process, 
the power curve turns around again at /i=1.832 [between 
'e' and T in Fig. [JJ (a)], and starts to decrease toward 
the first band [see 'g' in Fig. [JJ (b)]. Near the first band, 
the power curve turns around again, and moves back into 
the semi- infinite gap [see inset in Fig. [JJ (b)] . Near the 
turning point, the vortex spreads out again [see Fig. [5] 
(h,i)]. Shortly after turning back, the vortex exhibits an- 
other complex reshaping process evidenced by a twist in 
its power curve [see the inset in Fig. [JJ (b) near point 'i']. 
Afterwards, the power curve moves into the semi-infinite 
gap smoothly, and the vortex ultimately becomes four 
fundamental solitons with n/2 phase delay between each 



other. Fig. [8] (1) shows such a vortex when the propa- 
gation constant fi—1, which is far away from the Bloch 
band. The complexity of the power curve as well as the 
intricate reshaping process of vortices exhibited in this 
saturable-vortex family is very remarkable and is rarely 
seen in other wave systems. 

Now we consider off-site vortex solitons under focus- 
ing saturable nonlinearity. Again, many families of them 
are found. One such family is shown in Fig. [UJ (power 
curve) and Fig. [TO] (vortex profiles). We see that this 
power curve as well as its vortex reshaping process are 
simpler. Indeed, unlike the above on-site saturable case, 
the power curve here has no twists. It does have a turn- 
ing point near the Bloch band though, just like all other 
vortex families reported in this paper. At high powers, 
the off-site vortex has a ring profile [Fig. [JJJ (a)]. As its 
power decreases (with increasing [i values), it reshapes 
itself into a four-humped off-site vortex occupying four 
adjacent lattice sites [Fig. [TUJ (b, c)]. These vortices 
have been theoretically analyzed and experimentally ob- 
served in [U, [l^, Hfjj]. Comparing these vortices with 
those of off-site Kerr vortices in Fig. [31 we can see that 
this saturable-vortex family is the counterpart of the first 
off-site Kerr-vortex family. As the propagation constant 
/x moves from the lower branch to the upper one, this sat- 
urable vortex continuously reshapes itself, and eventually 
turns into four humps with 7r/2 phase delay between each 
other (as in the on-site case) . In all vortices of this off-site 
family, their phase fields have simple 2-7r-winding around 
the vortex center [see inset in Fig. [TUJ (a)], thus having 
charge one. 

When the saturable nonlinearity is of defocusing type 
(Eq < 0), infinite families of vortex solitons exist in the 
first band gap. In this case, we have found that saturable 
vortices show similar behaviors (such as power curves and 
vortex shapes) as defocusing Kerr vortices in Figs. [3JE1 
thus such results will not be displayed here. 



IV. SUMMARY AND DISCUSSION 

In this paper, we examined various families of charge- 
one vortex solitons in a 2D periodic media. For both 
Kerr and saturable nonlinearities, we found that there 
exist infinite families of on-site and off-site vortex soli- 
tons in the semi-infinite gap (for focusing nonlinearity) 
and the first gap (for defocusing nonlinearity). We fur- 
ther showed that all these vortex families do not bifurcate 
from edges of Bloch bands. Rather, before reaching band 
edges, they turn back and move into band gaps again. 
Within each vortex family, we revealed that topologi- 
cal structures of vortex solitons undergo drastic changes 
(such as evolving from a square configuration to a cross 
configuration) as the propagation constant varies. 

Regarding the linear stability of vortex solitons, it has 
been shown before that certain on-site and off-site charge- 
one vortex solitons in the semi-infinite gap (under focus- 
ing nonlinearity) are stable in deep lattice potentials or 
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(a) (b) (c) (d) (e) (f) 





FIG. 8: (Color online) On-site vortex profiles (|«|) corresponding to the marked points in Fig. [7] under focusing saturable 
nonlinearity. (1) shows the profile of the vortex at (far from the band edge). The inset in (a) is the typical phase structure 
of all these vortices. The background circles represent the lattice sites (with high / values), as in Fig. 1101 



with weak intersite couplings [H, [2(| H3, [H, • Certain 
on-site and off-site gap vortices under defocusing nonlin- 
earity are stable as well [1H H3| ■ Such vortices generally 
correspond to the lower branches of first-family vortices 
reported in this paper [see Figs. 3(a) and [UJb) for in- 
stance]. It will be interesting to investigate the linear 
stability of upper branches of first-family vortices as well 
as branches of higher-family vortices obtained in this pa- 
per. The stability results in [4Cj imply that some of the 
upper branches of vortices are linearly unstable [such as 
Fig. 2(d)]. In such cases, one can ask where on the vor- 
tex branches the instability starts to appear. Note that 
all vortices reported in this paper have charge one. It is 
known that vortices with more intensity humps can be 



stabilized by increasing their charges [I(J, 5l|. For m ~ 
stance, the vortex in Fig. 2(d) can be stabilized if its 
charge increases to three [4(| ■ Thus it would be interest- 
ing to study vortex families with higher charges as well 
as their stability properties. These questions will be left 
for future studies. 
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FIG. 9: The power diagram of an off-site vortex family in 
the semi-infinite gap under focusing saturable nonlinearity. 
Vortex profiles at the marked points are shown in Fig. 1101 
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FIG. 10: (Color online) Off-site vortex profiles (|u|) corre- 
sponding to the marked points in Fig. Sunder focusing sat- 
urable nonlinearity. The inset in (a) is the typical phase struc- 
ture of all these vortices. 
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